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Abstract 

A model of a carbon-oxygen (C-O) presupernova 
core with an initial mass 1.33Mq, an initial car- 
bon mass fraction — 0.27, and with an average 
mass growth-rate 5 • 1O _7 M0 yr _1 due to accretion 
in a binary system was evolved from initial density 
p c = 10 9 g cm~ 3 , and temperature T c = 2.05 • 10 8 K 
through convective core formation and its subsequent 
expansion to the carbon runaway at the center. The 
only thermonuclear reaction contained in the equa- 
tions of evolution and runaway was the carbon burn- 
ing reaction 12 C+ 12 C with an energy release corre- 
sponding to the full transition of carbon and oxygen 
(with the same rate as carbon) into 56 Ni. As a pa- 
rameter we take a c — a ratio of a mixing length 
to the size of the convective zone. In spite of the 
crude assumptions, we obtained a pattern of the run- 
away acceptable for the supernova theory with the 
strong dependence of its duration on a c . In the vari- 
ants with large enough values of a c — 4.0 • 10~ 3 and 
3.0 • 10 -3 the fuel combustion occurred from the very 
beginning as a prompt detonation. In the range of 
2.0 • 10~ 3 > a c > 3.0 • 10~ 4 the burning started as a 
deflagration with excitation of stellar pulsations with 
growing amplitude. Eventually, the detonation set 
in, which was activated near the surface layers of the 
presupernova (with m ~ 1.33A/0 ) and penetrated 
into the star down to the deflagration front. Excita- 
tion of model pulsations and formation of a detona- 



tion front are described in detail for the variant with 
a c = 1.0 • 10~ 3 . 

Keywords: supernovae and supernova remnants; 
plasma astrophysics; hydrodynamics and shock 
waves; detonation and deflagration. 

Introduction 

Evolution of a degenerate C-0 stellar core with a 
mass close to the Chandrasekhar limit (1.43M Q for 
a carbon white dwarf, see Bisnovatyi-Kogan, 1989) 
leads to the development of thermal instability and 
to the explosion. The critical density p c cr for the 
beginning of the explosion may vary from ~ 2 • 10° 
to ~ 10 10 g cm~ 3 . It increases with reduction of 
M — the core mass growth-rate which for a single 
asymptotic giant branch (AGB) star is defined by 
the Paczyhski — Uus relation (Paczyhski, 1970) 



M = 6 • 1Q- 7 (M/M e - 0.5)M Q yr" 



(1) 
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and can take values from 10~ 8 to 5 • 10 _7 M© yr _1 
in an accreting white dwarf (Iben, 1982; Hachisu et 
al., 1996). In the present paper we assume M = 
5 • 10~ 7 M Q yr" 1 which can be appropriate both for 
a single AGB-star and for a component of a binary 
system. 

The explosion calculations, some of which ac- 
counted for convection, were carried by various au- 
thors (Arnett, 1969; Ivanova et al., 1974; Nomoto 
et al., 1976; see also the review by Niemeyer and 
Woosley (1997)). Ivanova et al.(1974) found that 



1 



carbon was burning in the deflagration regime with 
excitation of pulsations, but convection at the super- 
nova stage was not considered, and the initial tem- 
perature profile was obtained not from evolutionary 
calculations, but from an estimate of a convection 
contribution. 

Woosley (1997) has carried out a series of 
one-dimensional explosion calculations for accreting 
white dwarfs (unfortunately, in this work the accre- 
tion rates are not specified) with critical central den- 
sities from 2 • 10 9 to 8.2 • 10 9 g cm" 3 and with nu- 
clear reactions for 442 isotopes, but he considered 
only the adiabatic convection at the presupernova 
stage, while at the beginning of the thermal run- 
away (when T c = 7 • 10 s K) he turned off the con- 
vection artificially (because further maintenance of 
the adiabatic temperature gradient in his model led 
to the prompt detonation), and defined the burning 
front propagation velocity according to Woosley and 
Timmes (1992) with a provision for a fractal dimen- 
sion of the front. 



Main equations of the convective 
hydro dynamic model 

Many authors (e.g. Iben, 1982) use the hydrostatic 
system of equations with adiabatic convection in the 
stellar core, when calculating the presupernova evo- 
lution (with time-scale ~ 10 4 yr) which is valid for 
this stage, but is too crude at the beginning of the 
thermal runaway. We carried our calculations both 
on the presupernova stage and on the stage of ex- 
plosive burning by the hydrodynamic scheme elab- 
orated in papers by Blinnikov and Rudzsky (1984), 
Blinnikov and Bartunov (1993), and by Blinnikov and 
Dunina-Barkovskaya (1993, 1994). In the latter two 
papers this hydrodynamic scheme was used for the 
calculation of white dwarf evolution. Here we have 
included non-adiabatic time-dependent convection in 
well-known mixing length approximation: 



dr/dt 
dv/dt 
dT/dt 



-Gm/r 2 -4irr 2 (dP/dm) , 
[e C c + £u + e g - 4wd(r 2 F conv )/dm 



(2) 
(3) 



dXc/dt 
du c /dt 



- 4nd(r 2 F Ilid )/dm\/(dE/dT) p , (4) 
= -X 2 r cc + (dX c /dt) conv , (5) 
= 2(v 2 - u 2 c )/l mix , (6) 



where Xq is the mass fraction of 12 C, ecc is the 
energy release during the carbon burning, e v stands 
for standard neutrino energy losses (as approximated 
by Schinder et al., 1987; see also Haft et al., 1994), 
e g is the energy liberation during the gravitational 
contraction, F ra( j is the radiative energy flux, F conv is 
the convective energy flux, u c is the time-dependent 
convective velocity (see Arnett 1969), and v c is the 
stationary convective velocity in the mixing length 
approximation which is given by the formula (see, 
e.g., Bisnovatyi-Kogan, 



v c = (g(dlnp/dr - {d\np/dr) s )) 1/2 • Z mix /2 ■ (7) 

The convective energy flux Fconv an d the value 
(dXc/dt) conv , which corresponds to the change of a 
carbon mass fraction Xq by convection, were calcu- 
lated by the formulas 



F r< 



= {{dT/dr) s - dT/dr)c p pu c l mi x/2 , (8) 



fdX c 



V dt 



r 2 dr 



r 2 l 



8X Q 
dr 



(9) 



The arrangement and the sizes of the convective zones 
are defined by the Schwarzschild criterion (according 
to the aforecited formulas (R) and (|h). On the presu- 
pernova stage we did not account for time-dependent 
convection and the formulas (||-§) contained v c in- 
stead of u c . A comparison of calculations with sta- 
tionary and time-dependent convection on the super- 
nova stage was performed in (Dunina-Barkovskaya, 
Imshennik, 2000). 

The mixing length in the i-th mass zone was de- 
termined by the relation 



mix 



= mm (ap(d\\nP\/dr) i 1 ,a c Ar c 



(10) 



where Ar c is the size of the convective zone contain- 
ing i-th zone, and ap in our calculations was taken 
equal to unity. 
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Thermodynamic functions for the electron- 
positron gas (the pressure P, the entropy S etc.) 
were calculated with help of Nadyozhin's asymptotic 
realtions described by Blinnikov et al. (1996), and 
for the ionic gas with screening — by Yakovlev and 
Shalybkov (1988). 

Initial and boundary conditions. 
Computations on the presuper- 
nova stage 

The initial state of the presupernova C-0 core should 
to be taken from reliable evolutionary calculations. 
We begin our computations with the central den- 
sity 10 9 g cm -3 which, according to Iben (1982), cor- 
responds to the central temperature about — 
2.05 • 10 8 K. Everywhere below we define, as a be- 
ginning of the runaway and, thereafter, as an end of 
the presupernova stage, the moment of time when the 
temperature in the central zone reaches the value of 
5 • 10 9 K. 

The initial carbon mass fraction Xq = 0.27 (con- 
stant throughout the C-0 core) was taken the same 
as in Iben's (1982) work. Assuming the above- 
mentioned central parameters for the C-0 core and 
its total mass 1.33M0, we have computed the adia- 
batic configuration in hydrostatic equilibrium which 
was used as an initial condition. It should be men- 
tioned that its difference from Iben's (1982) model 
seems insignificant. 

During the calculations of the evolution of this 
model with account for a constant mass growth-rate 
5 • 10~ 7 M Q yr" 1 (typical for accreting white dwarfs 
in binary systems), it is worth-while to set a non- 
zero outer boundary condition for the pressure. In 
our previous work (Dunina-Barkovskaya, Imshcnnik, 
2000), where the calculations led to relatively small 
changes of the total mass before the runaway, we in- 
cluded all the accreting mass in the outer boundary 
condition. 

This condition is easily derived in the approxima- 
tion of a thin (AR/R < 1) and light (AM /M < 1) 
envelope. Really, by integrating the hydrostatic equi- 
librium equation — dP/dm = Gm/Airr A over the en- 



velope thickness we obtain for the pressure Pb on the 
outer boundary of the zone with the Lagrangean co- 
ordinate M: 



Pb(M) = 



G_ f M+AM mdm 
4tt J m r 4 
GM 



G 



8ttP 4 



/ ovM+AM 



> AI 



4ttP 4 



AM, 



(11) 



where AM = Mt and R = R(t) during the hydro- 
static evolution, i.e. the outer pressure Pb from ( pl| ) 
increases in proportion to the time t which is mea- 
sured from the beginning of the C-0 core evolution 
calculations. 

In this case the full mass Mat + AM (where Mm 
is the mass of N Lagrangean zones contained in the 
hydrodynamic model) grew up to ~ 1.351AfQ, in ac- 
cordance with evolution time 4.25 ■ 10 4 yr (Dunina- 
Barkovskaya, Imshennik, 2000). 

In the present work, we have gradually, according 
to the accretion, added new Lagrangean zones to the 
model (cf. Woosley, 1997). The mass of the zones 
decreased geometrically from 7.83 • 10 _3 Af Q (151-th 
zone) down to 7.50 • 10~ 5 M Q (the last, 170-th zone). 
It was found in this case that the equation (|ll]) is 
quite inaccurate (evidently because AR/R is not very 
small). As a result, the evolution time on the presu- 
pernova stage increased nearly twice. For the vari- 
ants with non-zero a c the mass of the model reached 
1.3658M Q after t ~ 7.2 -10 4 yr, and the number of the 
Lagrangean zones became N = 170. After that we 
have added no Lagrangean zones, but simulated ac- 
cretion by increasing the outer pressure. In the run- 
away point (after t ~ 7.9 ■ 10 4 yr since the beginning 
of calculations) the outer pressure corresponded to 
the additional external mass from 2.316 -1O _3 M0 (for 
a c = 3.0-10- 4 ) to 2.342-10~ 3 M Q (for a c = 4.0-10~ 3 ). 
In the control variant with a c = the temperature 
in the central zone increased to 5 • 10 9 K already af- 
ter 6.789 • 10 4 years, therefore the model mass could 
grow only up to 1.3636M0, which corresponds to the 
number of zones N = 161. 

Let us consider the boundary condition ( pi] ) during 
the explosion for the variant with a c — 1.0 • 10 -3 
when AM = 2.338 • 10- 3 M Q , and R = 1.85 • 10 8 cm 
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(Fig. 1): 

PU1.33M )=P O 0^) . (12) 

where Pq = 5.72 • 10 22 din cm -2 . This value is to 
be compared with the central pressure of the C- 
O core which has the well-known lower estimate 
P c (0) > GM 2 /(87ri? 4 ). Thus, the ratio of the outer 
boundary pressure to the central one, Ph(M) / P c (0) 
< 2AM /M = 3.42 • 1CT 3 , is surely very small at the 
beginning of the runaway. In the subsequent runaway 
process it certainly changes and can increase by fac- 
tor of a few, but nevertheless remains small enough 
which can be seen from the numerical results. 

It is easy to estimate that the relative contribution 
of the inertial term into the derived boundary condi- 
tion (|lj) with characteristic parameters of C-0 core 
pulsations (4u p /t p ) • (fi^/(GM)) will be small: no 
more than a few percent. Thus we can neglect it in 
( |TT| ) even when the pulsations amplitude is maximal 
Ai? ~ Rq and r p ~ 5c by the end of deflagration 
burning (See Fig. 4). This boundary condition im- 
posed by the mass accretion onto the C-0 core sur- 
face is in close agreement with the outer boundary 
condition assumed in the previous works (Ivanova et 
al., 1974; Ivanova et al., 1977a) with the same depen- 
dence on an outer radius R : P\ } {M) oc i? -4 (p^). Let 
us remind that in the cited works the outer pressure 
also simulated the presence of the outer C-0 core en- 
velope. It was taken small enough, but about a fac- 
tor of 2 higher than the value from (|l2|). Thus, the 
effect of the outer boundary was more pronounced 
there than in the present work, but less significant 
than in (Dunina-Barkovskaya, Imshcnnik, 2000). It 
should be stressed that the non-zero outer pressure 
Pb(M) 7^ violates the energy conservation law dur- 
ing the explosion (see Fig. 4, 5) when the pulsations 
arise and the radius R of the C-0 core changes. 
Naturally, accounting for the accreting matter in the 
outer boundary condition (|ll]) and ( p"2| ) we imply its 
attachment to the evolving C-0 core without any fur- 
ther energy release, i.e. merely an adhesion. In fact, 
this substance has, generally speaking, another chem- 
ical composition, and can go through the thermonu- 
clear reactions with an energy liberation and even 



with partial ejection of the accreted material. We 
simply neglect those processes. On the other hand, 
the role of the outer boundary condition ( p"2| ) is very 
significant during the explosion, especially when the 
detonation wave forms near the surface (see below). 

The discussion of the results 

In the present work we have systematically inves- 
tigated the models numerically with different val- 
ues of the mixing length parameter ( J10| ) equal to 
a c = 4.0 ■ 10" 3 ; 3.0 ■ 10" 3 ; 2.0 • lO^ 3 ; 1.0 • 10~ 3 ; 
3.0 ■ 10~ 4 . From a physical standpoint, it is very diffi- 
cult to choose a certain value a c from this wide range 
in the used approximate theory of non-adiabatic con- 
vection. The energy release in the thermonuclear re- 
action of carbon burning 12 C+ 12 C -^ 24 Mg+7 was 
taken, first, by a rather crude formula obtained by 
Fowler and Hoyle (1965) with electron screening fac- 
tor by Salpeter (1954) (see also Salpeter and Van 
Horn 1969, and Arnett 1969), and, second, with a 
maximum possible energy liberation corresponding 
to the instantaneous combustion of all the carbon- 
oxygen mixture to 56 Ni. The energy release accepted 
here was taken from our early work (Ivanova et al. 
1974). It is convenient for the further comparison 
of our previous results with those obtained in the 
current work. Such a simplification of a thermonu- 
clear energy generation cannot affect essentially the 
obtained results due to the comparatively small sen- 
sitivity of the ignition conditions for the C-0 mixture 
to all pre-exponential values in the expression for the 
thermonuclear burning rate. 

As a result of the presupernova evolution calcu- 
lations we obtained that just prior to the runaway 
the density drops to a value which depends non- 
monotonically on the parameter a c and is maximal 
when a c = 4.0 • 10 -3 . In general, the initial critical 
central density (in units 10 9 g cm~ 3 ) p c cr 9 ~ (1.88 — 
2.03) turned out to be close to the one accepted ear- 
lier by Ivanova et al. (1974) and later by Ivanova 
et al. (1977a) (p c g = 2.33), though it is somewhat 
lower. Let us remind that in 1970s the most popu- 
lar viewpoint was that single intermediate mass stars 
evolve along the convergent Pazcyhski's (1970) track 
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Figure 1: The presupernova evolution: the depen- 
dence of the central temperature T c , the central den- 
sity pc, the convective core mass M conv , and the ra- 
dius Rn of the N-th (last) Lagrangcan zone upon 
time for a c = 1.0 • 10~ 3 after the convective core for- 
mation. 
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Figure 2: The temperature profiles at the begin- 
ning of the SN stage (when the central temperature 
reaches 5 • 10 9 K) for various values of the parameter 



after the formation of the carbon (or carbon-oxygen) 
core. The modern approach to the presupernova evo- 
lution involves the evolution in a close binary system 
with an accreting white dwarf. In this case, accord- 
ing to Yungclson (1998), all the stars in the initial 
main sequence mass range 2.5 < M/Mq < 10 may 
become carbon-oxygen white dwarfs in close binary 
systems. This range does not differ significantly from 
the similar range for single stars 3.5 < M/M Q < 8 
considered previously (Pazcyfiski, 1970). 

Some quantities characterizing the C-0 stellar core 
before the runaway are represented in Fig. 1 : the con- 
vective core mass M conv , the central temperature T c , 
the central density p c , and the radius of the C-0 core 
i?jv over the evolution time from 6 • 10 4 yr (the be- 
ginning of the convective core formation) to 7.87- 10 4 
yr for a c = 1.0 • 10~ 3 . Small jumps in the radius 
of the last (iV-th) Lagrangean zone shown on Fig. 2 



clearly indicate the procedure of adding of new La- 
grangean zones described above. The phase of the 
slow expansion before the runaway (see above) is al- 
most unnoticeable on this graph because its duration 
(~ 3 • 10 2 yr) is small in comparison to the overall 
evolution time (during which the star, as a whole, 
was contracting very slowly) and practically does not 
depend on a c . 

In Fig. 2, the initial profiles of temperature de- 
pending on the mass coordinate are represented at 
the moment of the beginning of explosion (see above). 
It can be seen that the temperature in the convective 
variants is appreciably higher than that in the vari- 
ant without convection. Even for a c — 3.0 • 10 -4 the 
temperature exceeds this value almost by a factor of 
2 in the central part of the convective core. Let us 
notice that in this moment there is an active combus- 
tion only in the first mass zone. On Fig. 2 we com- 
pare the initial temperature profiles with data from 
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Figure 3: The time dependence of the central temper- 
ature T c and of the full energy of the C-0 core E tot 
for a c = 4.0- 1CT 3 and 3.0 -10~ 3 (prompt detonation). 

Ivanova ct al. (1977b). In the latter work the explo- 
sion developed for the two (upper) T-profilcs. There- 
fore, it is not surprising that the runaway led to the 
total disruption of the star for all a c > 3.0 • 10~ 4 in 
the present calculation (see below). For the variants 
with a c = 2.0 • 10~ 3 and with a c = 3.0 • 10~ 4 which 
are the boundary values for the deflagration regime 
with pulsations (see below), the adiabatic tempera- 
ture profiles are also shown, with the entropy equal to 
the entropy of the 2nd mass zone for each variant. It 
can be seen that our profiles lie quite below the cor- 
responding adiabatic profiles, and they have another 
shape for m < O.1M , which of course can affect the 
development of the runaway without the regime of 
spontaneous burning (Blinnikov, Khokhlov, 1986). 

It is interesting that in the second mass zone (with 
m = 2.24 • 1O~ 4 M0) the smooth temperature pro- 
file for the maximum a c = 4.0 • 10~ 3 touches the 
third (lower) steep profile for which the explosion in 
(Ivanova et al., 1977b) failed. But the initial pro- 
file for the variant with the minimal a c = 3.0 • 10~ 4 
which in fact "hardly" exploded (see Fig. 5) already 
lies noticeably below the third profile from (Ivanova 



Figure 4: The time dependence of the central temper- 
ature T c and of the full energy of the C-0 core E tot 
for a c = 2.0 • 10~ 3 and a c = 1.0 • 10~ 3 (pulsational 
deflagration with delayed detonation) . 

et al., 1977b) for m < 7 ■ 10~ 4 Af Q (in the first 4 mass 
zones). We can conclude that there is a qualitative 
agreement between our present results and the calcu- 
lations in the cited work. 

The development of the runaway was different for 
a c > 3.0-10- 3 (Fig. 3) and for a c < 2.0-10~ 3 (Fig. 4). 
On the said graphs the time dependence of the cen- 
tral temperature T c and of the full energy of the C-0 
core -E^ot on time is shown. The latter eventually be- 
comes equal to the value ~ 1.5-10 51 erg corresponding 
to the full incineration into 56 Ni for all the variants 
regardless of the value a c . The value of E to t at the 
moment of the runaway beginning (—4.3 • 10 50 erg) is 
also practically the same for all this variants. 

The basic result of our numerical calculations is 
that the pulsational deflagration regime eventuates 
when the mixing length is small enough, namely 
a c < 2.0 • 10~ 3 . It was found that for a c > 3.0 • 10~ 3 
the prompt detonation occurs (see Fig. 3). Hereby 
we obtained a reasonably accurate threshold for the 
transition from the detonation to the pulsational de- 
flagration. 
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Figure 5: The time dependence of the central tem- 
perature T c and of the full energy of the C-0 core 
i?tot f° r a c = 3.0- 10~ 4 (pulsational deflagration with 
delayed detonation). 



On Fig. 4 there is 10 pulsations for the variant 
with a c = 2.0 • 10~ 3 and 14 pulsations for the vari- 
ant with a c = 1.0 • 10~ 3 . In both cases the burning 
ends with a powerful energy release in the last pul- 
sation. In this case the abrupt change of the burn- 
ing regime occurs. In the previous pulsations there 
was a deflagration regime, and in the last pulsation 
the all remaining fuel (about 90%) is burning in the 
regime of detonation. The value T c reaches its maxi- 
mal value (2.5 — 2.9) • 10 10 K a bit later, immediately 
after the collision of the detonation front propagat- 
ing inward from the surface and deflagration front, 
during the afterburning in the central zone). Thus, 
the hybrid regime of the thermonuclear explosion 
of a carbon-oxygen white dwarf takes place, which 
is most promising for the explosive nucleosynthesis 
(Niemeyer and Woosley, 1997). On the other hand, 
Fig. 3 demonstrates the burning without pulsations 
in an ordinary detonation regime starting from the 
center (for a c > 3.0 • 10~ 3 ). 

In Fig. 5, the time dependence of the central tem- 
perature and of the total energy is shown for the corn- 



Figure 6: The density profiles versus radius for a c = 
1.0 • 10~ 3 (pulsational deflagration). 

putation with a c = 3.0- 10~ 4 . Here the total number 
of pulsations is appreciably larger than on the previ- 
ous graph, but only 14 of them were accompanied 
by the propagation of deflagration into a neighboring 
mass zone, which can be well seen by very narrow 
specific maxima of the temperature T c . Let us notice 
that these narrow maxima correspond (due to the en- 
ergy generation) to the small jumps of the value -Etot- 
In the variant with a c = 1.0 ■ 10" 3 one to three mass 
zones were burning during each pulsation (except for 
5th and 7th), so almost all the pulsations had the 
said temperature peaks (Fig. 4). 

The deflagration regime is distinctly characterized 
by the profiles of density which are represented ver- 
sus the Eulerian radius r in Fig. 6 for a c = 1.0 • 10~ 3 . 
The density drops by a factor of 1.5 at the defla- 
gration front. Depending upon the pulsation phase 
the front either moves inward (the contraction phase) 
or outward (the expansion phase) during the overall 
propagation (by mass zones). The pulsations can be 
seen more clearly by the behaviour of the outer layers 
of the star, in the still unburnt matter. The ampli- 
tude of the outer radius pulsations grows in time and 
attains the value of ~ 1.5 (this can be seen well due 
to the logarithmic scale for density in Fig. 6). 
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Figure 7: The density and temperature profiles versus 
radius for a c — 3.0 • 10~ 3 (prompt detonation). 



The density profiles p(r) for stronger convection 
are given in Fig. 7 for a c = 3.0 • 10 ~ 3 . The profiles 
of the temperature T(r) are plotted on this graph in 
a small scale, they show neatly the position of the 
burning front as a jump of T. It can be seen that 
at the burning front the density of matter distinctly 
increases displaying the detonation burning regime. 
The supersonic character of detonation is also re- 
vealed by the quiescence of stellar layers outside the 
burning front. 

The convective processes are most neatly charac- 
terized by the profiles of the specific entropy which 
are shown in Fig. 8 for the case a c = 1.0 • 10~ 3 
versus the mass coordinate m (in logarithmic scale). 
The convective zones which are defined by the known 
Schwarzschild criterion, i.e. correspond to the regions 
with the negative entropy gradient, are marked on 
these profiles for some typical moments of time. This 
figure shows that the primary convective core during 
the runaway is being separated into some convective 
zones, the largest of which includes one of two mass 
zones behind the burning front and the rest of the 
convective core before the front. The separation pro- 
cess is of nonstationary character and indicates the 
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Figure 8: The entropy profiles versus mass for a c = 
1.0 • 10~ 3 (before the detonation). Convective zones 
are shown by labels and numbers, e.g. 2-1 means 
the first convective zone for the second moment 
(t=0.6222 s). 

importance of the time-dependent convection. 

The entropy profiles are given in Fig. 8 before the 
formation of the detonation wave propagating to- 
wards the center (see Fig. 4). After the formation 
inception of such a wave, naturally, in a very short 
time <~ 0.1 s, the entropy behind the detonation front 
rises to the values which are typical also for the defla- 
gration, (2 — 3) • 10 8 crg g _1 K _1 ), i.e. by factor of 4-6 
more than an initial entropy (~ 0.5 • 10 8 erg g _1 K _1 ), 
with a positive gradient. Only in outer layers (with 
m > 1.36-Mq) the negative entropy gradient persists, 
and, therefore, the convection occurs, partly prevent- 
ing the detonation process described here. 

The important role of the convection intensity 
characterized by a parameter a c attracts our special 
attention. It is clear that the deflagration burning 
regime with pulsations (of course, the deflagration is 
not necessarily accompanied by pulsations, but for 
the latter the subsonic deflagration is the necessary 
condition) takes place in the whole range of the a c 
values, a c m i n <; a c <; 2.0 • 10~ 3 . It is essential that 
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oi c min is above zero, because for a c = in our calcu- 
lation there was no hydrodynamic explosion at all, al- 
though the first mass zone was burning. We note that 
for higher critical densities p c cr which is the case for 
lower accretion rate (see above) , the value of a c m i n 
can drop practically down to zero (Zmitrenko et al., 
1978). However, for high densities p c cr the account 
for neutrino energy losses together with the neutron- 
ization kinetics is required (Ivanova et al., 1977a, c). 
But in the present case (p c cr ~ 2 • 10 9 g cm~ 3 ) the 
role of these complicated processes is not important, 
which was shown by direct comparison of the calcu- 
lations of Ivanova et al. (1974) with those of Ivanova 
et al. (1977b) at the close value of central density 
p c = 2.33 • 10 9 g cm~ 3 . Thus, the account for con- 
vection processes can prevent the pulsations develop- 
ment during the deflagration burning when the pa- 
rameter a c exceed some critical value a c cr ; t . In the 
present work it is found that a c cr ; t > 3.0 • 1CP 3 . 

A qualitative analysis of the 
physical conditions of the de- 
layed detonation 

In essence, the emergence of the detonation, after 
a long period of the deflagration burning accompa- 
nied by the development of the global pulsations of 
the C-0 core, can be termed a " delayed detonation" 
(Niemeyer and Woosley, 1997): from the moment 
of the central ignition a long time (tens of seconds) 
passes, i.e. many thermodynamic time-scales of the 
order of the pulsation period, 2-3 s. An analysis of the 
results has shown that the detonation burning front 
emerges in a mass coordinate about 1.33M Q , more 
specifiacally, in the 150th mass zone when the total 
number of zones equals 170. This coordinate is virtu- 
ally independent of the parameter a c . In our model in 
the 150th zone the mass grid changed from uniform 
to geometrically decreasing. In the work (Dunina- 
Barkovskaya, Imshennik, 2000) the 127th mass zone 
in which the detonation begins had the same feature, 
so we have to investigate the influence of a mass grid 
nonuniformity on a detonation initiation in our fu- 
ture work. In the present section we shall try to give 




1.29 1.31 1.33 1.35 1.37 



m, M Q 

Figure 9: The entropy profiles versus mass in the 
outer layers of the C-0 core for a c — 1.0 TO -3 during 
the delayed detonation (up to the beginning of last 
pulsation with the detonation initiation). 

a physical justification of the entropy growth that 
leads to the detonation. 

After some first pulsations, i.e. long before the 
generation of the burning front, the entropy of the 
outer mass zones begins to increase (Fig. 9). At the 
moment t — 45.098 s (the beginning of the last pul- 
sation) two local entropy maximums exist, the out- 
ermost one (to ~ 1.36M Q ) and the innermost one 
(to ~ 1.33M Q ). Though the entropy in the inner- 
most maximum is less than in the outermost one, 
the temperature in the innermost maximum during 
the last pulsation is higher (due to the higher den- 
sity). Just at this maximum the detonation burning 
front is born (see Fig. 11). We can see the relatively 
high growth of entropy between the close moments 
t = 2.577 s and t = 2.850 s, between the maximum 
and the minimum of the second (!) pulsation (see 
Fig. 4), which is less by amplitude than the subse- 
quent pulsations. 

It can be seen on the graphs for density (Fig. 10) 
that the moments close to the maxima of pulsations 
are accompanied by a subsequent emergence of the 
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m , M Q 

Figure 10: The density profiles versus mass in the 
outer layers of the C-0 core for a c = 1.0 ■ 10~ 3 (pul- 
sational deflagration up to the detonation initiation) . 



compression waves on the stellar surface. Follow- 
ing Ivanova et al. (1982), one can estimate the en- 
tropy growth AS* treating these compression waves as 
weak shocks with a pressure jump AP. The one gets 
the well-known entropy jump (Landau and Lifshitz, 
1954): 

1 fd 2 V s 



AS 



12Ti \dP? 



{Apy 



(13) 



This estimate is easy to apply in the case under 
consideration because we can account only for the 
pressure of the completely degenerate ultrarelativis- 
tic electron gas P\ c = Kp\^ , with K = 4.90 • 



-1/3 

abatic derivative in 



10 14 cm 3 g i/3 s a , in calculation of the second adi- 



11311 , whence it follows that 



fd 2 v\ _ f d 2 v 



21 



d 2 V 



dP 2 e 16 K 2 p n / 3 ' 



(14) 

Such an approximation is possible because in the zone 
under consideration with m = 1.33Mq we can take 
as an estimate for the values p\ and T± in front of 
the involved weak shock their values at the moment 
t = 2.577 s, namely, p x = 1.073 • 10 7 g cm" 3 , and 
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Figure 11: The temperature profiles versus radius 
during the delayed detonation for a c = 1.0 • 10 -3 . 

Ti = 6.64- 10 7 K. After a simple rearrangement we get 
the final expression for the specific entropy growth in 
(TL1) (S is in units 10 8 erg g _1 K-l, see Fig. 9): 



ASs = 0.536 • 10 



1.78 



1/3 
6 Pi 

Ti 

4/3 



El 
Pi 



- 1 



L± _1 



(15) 



In the expression (|T^) it is necessary to specify 
also the amplitude of the weak shock wave. From 
Fig. 10 we can get pi = 1.35 • 10 7 g cm -3 (for the 
moment t = 2.850 s). Then we finally get from ([l5|): 
ASs — 0.085, whereas in Fig. 9 the entropy growth 
between these two moments is AS$ ~ 0.06. Thereby 
we can reasonably explain the entropy growth ob- 
tained in calculations as a result of a weak shock dis- 
sipation which occurs in the region of a rather abrupt 
density decrease. Closer to the surface the entropy 
growth subsides, according to ([l5|), because of the 
density lowering towards the star edge, and inward 
it subsides too, probably due to the wave amplitude 
reduction (see also Fig. 10). Nevertheless, the local 
entropy maximum (Fig. 9) does not disappear with 
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time, but persists and increases furthermore by a fac- 
tor of ~ 2.7 before the moment of the detonation gen- 
eration. It reaches the value of 1.6 • 10 8 erg g _1 K— 1. 
Let us notice that unfortunately this cannot be seen 
on Fig. 8, because the outer stellar layer are shown 
there in the logarithmic scale by mass. 

In Fig. 11, we represent the temperature profiles 
from the moment of the burning front initiation 
(t = 45.098s) to its "collision" with the deflagration 
front at a radius r ~ 5 • 10 7 cm (which corresponds 
to the mass from the center 0.14M Q ) at the moment 
t = 45.212 s. On this graph the beginning of the 
expansion of matter behind the detonation wave can 
be well seen: the outer radius rises from 2.3 • 10 s cm 
to 3.5 • 10 s cm at the temperature ~ 10 9 K. When 
the detonation wave propagates in the direction of 
increasing density, the temperature rises to 10 10 K, 
partly due to the heat capacity decline. The com- 
pression on the front is typically modest: only by a 
factor of 1.5 (see Fig. 7). The pressure jump is not 
high too: about a factor of 2 till the end. It is re- 
markable that the generation of the detonation front 
really occurs at the innermost entropy maximum (see 
Fig. 9), which corresponds in Fig. 11 to the innermost 
temperature maximum at the moment t — 45.098 s. 

A fair question arises: whether such a detona- 
tion front is stable with respect to galloping insta- 
bility considered recently by Imshennik et al. (1999). 
This question requires an additional analysis, but it 
is conceivable that the stability here is more natu- 
ral, because the density in front of the detonation 
wave increases which counteracts the decoupling of 
the shock wave front from the burning zone. There 
is one more argument in favor of the stability against 
the multidimensional disturbances, namely, the large 
enough width of the burning zone which can be com- 
pared with the C-0 core radius at the densities 
p ~ 10 7 gcuT 3 (Imshennik, Khokhlov, 1984). This 
implies also that the burning of the C-0 mixture as 
a matter of fact does not fuse the elements to "iron" 
("Fe") in the converging detonation wave, but is lim- 
ited, for example, by such nuclides as Si etc. The 
problems of nucleosynthesis in the obtained regime 
of delayed detonation should be investigated in other 
works. 



Conclusion 

In the current study, first of all, we obtained the 
initial conditions, justified by evolution calculations, 
for the thermonuclear runaway in the carbon-oxygen 
stellar core with a mass close to the Chandrasekhar 
limit due to the accretion prescribed by a constant 
value 5 • 1O~ 7 M yr _1 , which docs not contradict 
to the current evolutionary studies of close binary 
systems. In our simulations the runaway begins at 
the moment defined by the calculation itself (and the 
equations do not change, except for the account for 
the time-dependent convection). Thus, we do not 
have to solve a complicated problem of the choice 
of the initial conditions, which largely affected the 
inception of the runaway. In our previous work 
(Ivanova et al., 1974, 1977a, b,c) the initial conditions 
(i. e. the initial temperature profiles) were defined, 
strictly speaking, without an adequate justification. 

The present computations include a well-known 
approximate model of a nonadiabatic convection, 
namely, the mixing length theory with the only ar- 
bitrary parameter a c (a ratio of the mixing length 
Z m i x to the radial extent of the convective zone Ar c ). 
The value of this parameter could be, in principle, 
extracted from the multidimensional approaches de- 
veloped currently (see Lisewski et al., 2000), i. e. the 
arbitrariness of the approximate theory could be re- 
moved. As a result of our computations it was found 
that the delayed detonation regime takes place under 
quite a wide range of values 3 • 10~ 4 £ a c £ 2 • 10~ 3 . 
The lower boundary of this interval is defined as an 
upper estimate. At higher values of a parameter 
a c > 3 • 10~ 3 in our calculations the ordinary detona- 
tion propagating from the center of the stellar C-0 
core was obtained, which is unlikely to take place be- 
cause of the instability of the burning front etc. It 
could be interesting to find out whether the specified 
interval intersects with the region of effective mix- 
ing length parameters justified by a multidimensional 
theory of turbulence. But already with the paramet- 
ric simulations of the flame propagation velocity we 
can impose some restrictions on the parameter a c . 
For example, we can notice that for a c = 1.0 • 10~ 3 
during the burning of a few first mass zones this ve- 
locity was close to the laminar front velocity obtained 
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by Woosley and Timmes (1992). 

In the present work, due to an application of 
the new hydrodynamic code with a variable La- 
grangean (difference) grid and with account for an 
outer boundary pressure conditioned by accretion, 
we have obtained a scenario of the development of 
a delayed detonation propagating from the surface 
towards the center of the star. The detonation is, in 
all probability, stable with respect to the galloping 
instability, and is conditioned by a preceding stage 
of the pulsational deflagration. This may prove to be 
important because of the recently emerged scepticism 
(Lisewski et al., 2000) in regard to the scenario of the 
deflagration to detonation transition due to the de- 
struction of the laminar burning front by turbulent 
eddies published earlier (Khokhlov et al., 1997). 

This work has been supported in part by the ISTC 
grant # 0370 and by the RFBR grant # 00-02-17230. 
The authors are grateful to W.Hillcbrandt for the 
helpful discussion. 
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